arXiv: 1506.0353Ivl [quant-ph] llJun2015 


Casimir—Polder force between anisotropic nanoparticles and gently curved surfaces 


Giuseppe Bimonte 

Dipartimento di Fisica, Universitd di Napoli Federico II, 

Complesso Universitario di Monte S. Angelo, Via Cintia, 1-80126 Napoli, Italy and 
INFN Sezione di Napoli, 1-80126 Napoli, Italy 

Thorsten Emig 

Lahoratoire de Physique Theorique et ModNes Statistiques, CNRS UMR 8626, 

Bat. 100, Universite Paris-Sud, 91405 Orsay cedex, France 
Massachusetts Institute of Technology, MultiScale Materials Science for Energy and Environment, 
Joint MIT-CNRS Laboratory (UMI 3466), Cambridge, Massachusetts 02139, USA and 
Massachusetts Institute of Technology, Department of Physics, Cambridge, Massachusetts 02139, USA 

Mehran Kardar 

Massachusetts Institute of Technology, Department of Physics, Cambridge, Massachusetts 02139, USA 

(Dated: June 12, 2015) 

The Casimir-Polder interaction between an anisotropic particle and a surface is orientation de¬ 
pendent. We study novel orientational effects that arise due to curvature of the surface for distances 
much smaller than the radii of curvature by employing a derivative expansion. For nanoparticles 
we derive a general short distance expansion of the interaction potential in terms of their dipolar 
polarizabilities. Explicit results are presented for nano-spheroids made of Si02 and gold, both at 
zero and at finite temperatures. The preferred orientation of the particle is strongly dependent on 
curvature, temperature, as well as material properties. 
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I. INTRODUCTION 

The interaction of small particles with surfaces is im¬ 
portant to a plethora of phenomena in physics, chemistry 
and biology. While the cause of the interaction can differ, 
in many situations the particles are neutral (source free) 
and their interaction with the surface is due to an embed¬ 
ding fluctuating medium or field. There is considerable 
interest in investigating how the interaction is affected 
by the geometrical shape of the surface, and several ex¬ 
periments have probed dispersion forces between 

particles and micro-structured surfaces. More elaborate 
examples for this type of interaction include quantum 
frictional forces acting on particles moving along a sur¬ 
face [5], the heat transfer between nanoparticles and 
curved or rough surfaces [6], and critical Casimir forces 
in colloidal systems and superfluid helium 13 El. 

Here we consider forces induced by quantum (and 
thermal) fluctuations of the electromagnetic (EM) field, 
known as van der Waals or Casimir-Polder interactions. 
The forces between a particle (atom) and a flat sur¬ 
faces have been extensively studied [9] ; for recent reviews 
see However, roughness and curvature which are 

ubiquitous features of many surfaces modify fluctuation- 
induced forces. Computing such interactions is compli¬ 
cated by their characteristic non-additivity which leads 
to interesting effects for anisotropic particles [13]. Eor 
good conductors, the force between spheroids scales not 
with the product of their actual volumes but with the 
product of the volumes of the enclosing spheres [T4| . The 


classic result of Balian and Duplantier m and more re¬ 
cently developed scattering techniques nnun] have been 
most successfully applied at distances that are large com¬ 
pared to the radii of curvature of the surface, and for a 
few specific surface shapes. A perturbative approach is 
presented in [18], where surfaces with smooth corruga¬ 
tions of small amplitude, were studied. The validity of 
the latter is limited to particle-surface separations much 
larger than the corrugation amplitude. 

However, the regime most relevant to experiments is 
at short distances (compared to the radii of curvature 
of the surface). Analytical results are known only for 
specific geometries, like for a perfectly conducting cylin¬ 
der and an atom m- A commonly used method in this 
regime is the proximity force approximation (PEA) [12], 
based on integrating the force to a fiat plate over vary¬ 
ing separations. This approximation clearly fails for 
anisotropic particles whose preferred orientation depends 
on the shape of the nearby surface. Here, we employ 
a systematic approach that becomes exact in the limit 
of small particle-surface separations. It is based on an 
expansion of the interaction potential in derivatives of 
the surface profile, and hence applies to general, curved 
surfaces. An analogous expansion has been used re¬ 
cently [20H22] to study the Casimir interaction between 
two non-planar surfaces. It has also been applied to other 
problems involving short range interactions between sur¬ 
faces, like radiative heat transfer m, and stray electro¬ 
static forces between conductors [24] . 

The paper is organized as follows: In Sec. H we present 
the derivative expansion for the general case of a particle 
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with electric and magnetic dipolar polarizabilities in front 
of a dielectric curved surface, and we specialize the results 
to the perfectly reflecting limit. In Sec. Ill we compute 
explicitly the orientation dependence of the interaction 
for spheroids made of Si 02 and gold, both at zero and at 
finite temperatures. Section IV summarizes our results 
and provides an outlook. 

II. DERIVATIVE EXPANSION OF THE 
CASIMIR-POLDER POTENTIAL 

Consider a nanoparticle near a dielectric surface S. We 
assume that the particle is small enough (compared to 
the scale of its separation d to the surface) to be be con¬ 
sidered as point-like, with its response to the electromag¬ 
netic fields fully described by the electric and magnetic 
dipolar polarizability tensors and respec¬ 

tively. Let us denote by the plane through the particle 
which is orthogonal to the distance vector (which we take 
to be the z axis) connecting the particle to the point P 
of S closest to the particle. We assume that the sur¬ 
face S is characterized by a smooth profile ^ = ^(^, 
where x = {x^y) is the vector spanning Ei (see Fig. [^. 
In what follows Greek indices /i, z^, • • • label all coordi¬ 
nates (x, y, z), while latin indices i, j, /c, • • • refer to (x, y) 
coordinates in the plane Ei. Throughout we adopt the 
convention that repeated indices are summed over. 

The exact Casimir-Polder potential at finite tempera¬ 
ture T is given by the formuia nsun] 

oo 

U = -A:BT^'Tr[T(®)UT(^)U](K„) . (1) 

n=0 

Here and denote, respectiveiy, the scattering 
T-operators of the piate S and the particie, evaiuated at 
the Matsubara wave numbers Kn = 27rnkBT/{hc)^ and 
the primed sum indicates that the n = 0 term carries 
weight 1/2. In a piane-wave basis |k, Q) |25] where k 
is the in-plane wave-vector, and Q = E^M labels re¬ 
spectively electric (transverse magnetic) and magnetic 
(transverse electric) modes, the translation operator U 
in Eq. Q is d iagonal with matrix elements where 
q = = g(/c), k = |k|. The matrix elements of 

the particle T-operator in dipole approximation are 

rQ3(k,k') = 

+ ’ ( 2 ) 

where q' = q{k'), e^^(k) = zxk, e^^(k) = — l//i:^(i/czzb 
q'k), k = k//c and we set 


The T-operator of an arbitrary curved plate is not 
known in closed form, and its computation is in general 
quite challenging, even numerically. In Ref. [26], how¬ 
ever, the leading curvature corrections to the potential 
were computed for an atom in front of a smoothly curved 



FIG. 1: Parametrization the configuration of a nano-spheroid 
near a gently curved surface: Solid curves on the surface S 
indicate the principal directions at P with local radii of cur¬ 
vature Ri and R 2 . Positive (negative) Rj correspond to a 
surface that curves away from (towards) the particle. 


surface, in the experimentally relevant limit of small sep¬ 
arations. The key idea is that as the Gasimir-Polder in¬ 
teraction falls off rapidly with separation, it is reasonably 
expected that the potential U is dominated by a small 
neighborhood of the point P of S which is closest to the 
particle. This physically plausible idea suggests that for 
small separations d, the potential U can be expanded 
as a series in an increasing number of derivatives of the 
height profile id, evaluated at the particle’s position. Up 
to fourth order, and assuming that the surface is homo¬ 
geneous and isotropic, the most general expression which 
is invariant under rotations of the (x, y) coordinates, and 
that involves at most four derivatives of H (but no first 
derivatives since Vid(O) = 0) can be expressed (up to 
0{d~^)) as 
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where + ayy^ and it is understood that all 

derivatives of H{x.) are evaluated at the particle’s posi¬ 
tion, i.e., for X = 0 . The coefficients are dimension¬ 
less functions of = ‘I'^nksT/{hc)^ and of any other 
dimensionless ratio of frequencies characterizing the ma¬ 
terial of the surface. The derivative expansion in Eq. @ 
can be formally obtained by a re-summation of the per¬ 
turbative series for the potential for small in-plane mo¬ 
menta k [26]. We note that there are additional terms 


involving four derivatives of H which, however, yield con¬ 
tributions ^ 1/d (as do terms involving five derivatives 
of H) and are hence neglected. 

A geometrical interpretation of Eq. (§ is obtained 
when the x and y axis are chosen to coincide with the 
principal directions of curvature of S at P. Then the 
expansion of is = d + x‘^/{2Ri) + y‘^/{ 2 R 2 ) + • • •, 
where Ri and R 2 are the radii of curvature at P. In this 
coordinate system, the derivative expansion of U reads 


U = - 


ksT 


{ (2) 
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As demonstrated in Ref. [26], the coefficients in 
Eq. © can be extracted from the perturbative series of 
the potential U. To second order in the deformation 
h(x) = P(x) — d, this involves an expansion of the T- 
operator of the surface S to the same order. The latter 
expansion was obtained in Ref. [27] for a dielectric ma¬ 
terial described by a frequency dependent permittivity 
e(co’). It reads 


Tqq, (k, k') = (27r)2(5(2) (k - k') Sqq, {icK^, k) 

"-2BQQ,(k,kOMk-kO 
(i? 2 )QQ'(k,k';k")/i(k-k")ft(k"-k') + 


+ 

"d^k" 


( 5 ) 




k and k' (as it is apparent from Eq. 0 ) that cannot be 
performed analytically for a dielectric plate. In the fol¬ 
lowing, we shall consider a perfect conductor, in which 
case the integrals can be carried out analytically. In this 
case, the matrix PQQ/(k, k') takes the simple form 


i?(k,k') 


k-k'/t^+/c/e' 

^ I 

• (k X k') 


^z.(kxk')^ 
-k • k' y ’ 


(6) 


where the matrix entries Q, Q' correspond to P, M re¬ 
spectively. Also, the matrix (P 2 )qq'( k, k'; k") is simply 
related to B by 


where (ic/^n, k) denote the familiar Eresnel reflection 

coeffcients of a flat surface, and h(k) is the Eourier trans- 
fromed deformation. Explicit expressions for the kernels 
PQQ/(k, k') and {B2)qq0^', k'; k") are given in Ref. [27]. 
Computing the coefficients f3p^^ involves an integral over 


(P 2 )(k, k'; k") = 2 g"P(k, k")a 3 P(k", k') , (7) 

where (73 = diag(l,—1). The coefficients are now 
functions of ^ only, and we list them in Tables [I| [II] for 
electric and magnetic dipole polarizabilities, respectively. 
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4 

1 

^(3 + 6 C + 15e + 22 f + 2^4 - 4f) 



2 

- 910(15 + 542^ + 259^^ - 546^’' - 14^^ + 28^®) 

- 2 e^(i-^ + Mo) 


3 

^{15 + 30^ - + 70C® + 2 C^ - 4f) 

§4(18-a 


4 

^(45 + 218^ - 59^^ + 146^’' + 14^^ - 28^®) 

lo(40-7e^) 


5 

^(9 + 18e - 27^^ + 50^’' - 2 ^^ + 4^®) 

(1 + S) 


TABLE I: The coefficients for the electric dipole contribution are obtained by multiplying the third column by e and 
adding the fourth column times Ei(2^) = — dtexp(—t)/t. 
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4 

- 4 (3 + 6^ + 33^2 - 70^3 - 2^^ + 4^®) 

Ml ('i. 

2 18 J 


5 

- jlg(15 + 158^ + 121^2 - 214^3 + 14^^ - 28^®) 

2 3 150 J 


TABLE II: The coefficients for the magnetic dipole contribution, using the same notation as in Tab. ffi 


III. ORIENTATION DEPENDENCE 


In this section we investigate the shape and orientation 
dependence of the Casimir-Polder force using Eq. Q. 
Before we consider a curved surface, it is interesting to 
stress that for a perfectly reflecting planar surface there is 
no orientation dependence at zero temperature for dipo¬ 
lar particles with frequency independent polarizabilities 
as realized, e.g., in the perfectly conducting limit. This 
follows directly from the fact that the ^ integrals of the 
two coefficients and /3p|2 equal so that the po¬ 
tential is proportional to the rotationally invariant trace 
of [14]. Coming back to a curved surface, we 

assume for simplicity that its height profile H is invari¬ 
ant under independent reflections in the x and y direc¬ 
tions. This symmetry of the surface ensures that the term 
proportional to /3p^ in Eq. Q is absent. Moreover, we 
assume that the particle has one axis of rotational sym¬ 
metry. In a new orthogonal basis (1,2,3) oriented such 


that the third axis coincides with the particle’s symme¬ 
try axis C, the polarizability tensors are diagonal with 
= diag(a^/2,a^/2,d|’3). 

The polarizability tensors for an arbitrary orientation 
are then obtained as a = lZ~^alZ^ where IZ is the matrix 
that rotates the principal axis of the particle to the basis 
composed of the principal directions of the surface S', i.e. 
7^(1, 2, 3) ^ {x^y^z). The orientation of the particle is 
conveniently parametrized by the polar angles (^, 0) of its 
symmetry axis C, where 0 is the angle formed by C and 
the z axis, and 0 is the angle between the (z, C) plane and 
the (x^z) plane (see Fig. [^. The polarizability tensors 
in the two coordinate systems are related by 

~ ^ ^ ^ ~ cos(26>)] , (8) 

o^zz = ^ [«± + 2 dfa + cos(26»)] , (9) 

(xL -0‘vy = ^ cos(2(/)) sin^ 0 , (10) 
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where we defined = 2 . Since the x and y 

axis are chosen to coincide with the principal directions 


of S at P, Eqs. (8pQ) together with Eq. Q yield the 
potential in the simple form 


U = - 


ksTV 

(P 


A{d) + B{d) cos( 2 ^) + C{d) ( ] cos(20) sin^( 6 >) 

Pi P2 


( 11 ) 


where A{d), C{d) are in general functions of tem¬ 

perature and the ratios d/P^, but do not depend on the 
angles 0 and 0, and V is the volume of the particle. Be¬ 
fore turning to detailed computations, we briefly discuss 
the qualitative features of the potential. It is obvious 
that the coefficients B and C must vanish for a spherical 
particle. Eor a non-spherical particle in front of a planar 
surface (Pi = P 2 ^ oc) the potential U in Eq. is 
invariant under rotations about the z-axis, as expected. 
The coefficient P, however, is in general different from 
zero, and hence even for a planar surface the potential 
depends on the polar angle 0 (except, as discussed above, 
for a perfectly reflecting surface and for frequency inde¬ 
pendent polarizabilities). In order to have a non-trivial 
dependence of U on the azimuthal angle 0, it is neces¬ 
sary to break the rotational symmetry about the 2 ;-axis. 
Clearly, this happens when the surface has different radii 
of curvature Pi 7 ^ P 2 at P, as evidenced by the third 
term between the brackets of Eq. 0 - Eor Pi 7 ^ P 2 , 
it is easy to verify that in general the potential U has a 
unique minimum, corresponding to an orientation of the 
particle along one of the (x, y, z) axes. More precisely, 
with D = CidjRx — djRP^^ the stable orientation of the 
particle’s symmetry axis is along the 



FIG. 2: Stable orientations of the particle’s symmetry axis 
as function of the coefficients B and D = C{d/Ri — d/R 2 ) in 
Eq. In). 


(i) X axis if P > max{0, 2P}, 

(ii) y axis if P < min{0, —2P}, 

(iii) z axis otherwise. 


in Eq. Q the expressions of the coefficients for a 
perfect conductor, which are listed in Tables I and II. 


The stable orientations are summarized in the diagram 
of Eig. To grasp more easily the different orientations 
of the particle relative to the curved surface, we show in 
Eig. the typical surface shapes for positive and negative 
radii of curvature, along with the coordinate frames for 
the position and orientation of the particle. 

Below, we numerically compute the potential between 
a gold surface and a spheroidal particle, made either 
of gold or of vitreous Si 02 . Eor particle-surface sep¬ 
arations d larger than the plasma wavelength of gold, 
Xp = 27rcIujp 120 nm, and smaller than a few micron, 
as we shall consider, the penetration depth of the elec¬ 
tromagnetic fields in gold contributing to the Casimir- 
Polder potential is (igoid < 20nm, and therefore it is al¬ 
ways much smaller than the separation d. In this range 
of separations, the gold surface S can thus be considered 
as perfectly reflecting, and it is therefore justified to use 


A. Si02 particle 

Eor a dielectric ellipsoid with electric permittivity e 
(and magnetic permeability /i = 1 ), the polarizability 
tensor is diagonal with respect to its principal axes, 
with elements (for y G {1, 2, 3}) 

~F ^ e — 1 

where V = 47 rrir 2 r 3/3 is the ellipsoid’s volume. In the 
case of spheroids, for which vi = r 2 = R and rs = T/ 2 , 
the so-called depolarizing factors can be expressed in 
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FIG. 3: Typical surface shapes for limiting cases correspond¬ 
ing to combinations of i?i = ±1, i ?2 = ±3 and R 2 ^ oo (ar¬ 
bitrary units). The coordinate frames indicate the position 
and orientation of the particle. 


terms of elementary functions, 

1 - ns 

m = n2 = —T— 




( 13 ) 


where the eccentricity e = jL‘^ is real for a 

prolate spheroid {L > 2R) and imaginary for an oblate 
spheroid (L < 2R). For a prolate spheroid 0 < ns < 1/3, 
while for an oblate spheroid 1/3 < ns < 1, the value 
ns = 1/3 corresponding to a sphere. For the dynamic 
permittivity e(ic(;) along the imaginary frequency axis, 
we use the simple two-oscillator model 


5(ia;) = 1 + 


CuV^UV , 


+ 


(14) 


‘"UV 


''IR 


with the parameters Cuv = 1.098, Cir = 1.703, cjuv = 
2.033 X 10^^ rad/s, and cjjr = 1.88 x 10^^ rad/s, which 
were obtained by a fit to optical data for Si 02 [28] . 

We observed earlier that the potential U is minimized 
when the particle’s axis points in the direction of one of 
the coordinate axes. First we determine the preferred 
orientations at zero temperature. In Fig. we show 
the stability diagram for a Si 02 oblate spheroid with 
ns = 0.7 (“pancake”) and Ri = 1000 /im. An analogous 



FIG. 4: Stability diagram for a Si02 oblate spheroid (“pan¬ 
cake”) with ns = 0.7 and Ri — lOOO/xm at T = 0 K. 



z-axis 



0.00 

d/R, 


FIG. 5: Stability diagram for a Si02 prolate spheroid (“nee¬ 
dle”) with ns = 0.2 and Ri = 100 /xm at T = 0 K. 


diagram for a prolate spheroid with ns = 0.2 (“needle”) 
and Ri = 100 /xm is show in Fig. From the diagrams 
it can be observed that the signs of surface curvature 
have an important effect on the preferred orientation of 
the particle. We note that the diagram depends on the 
choice of Ri since it is compared to the material depen¬ 
dent length scales that are set by the characteristic fre¬ 
quencies of Si 02 . 

Thermal fluctuations have a strong impact on the sta¬ 
ble orientation of the particle. For room temperature, 
T = 300K, the stability diagram for a “pancake” is shown 
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range of the stability plots shown here. 





FIG. 6: Stability diagram for a Si02 oblate spheroid (“pan¬ 
cake”) with ns = 0.7 and Ri = 1000/xm at T = 300 K. FIG. 8: Stability diagram for a gold prolate spheroid (“nee¬ 

dle”) with ns = 0.2 at T = 0 K. 



d/Ri 


B. Gold particle 


Next we consider a gold spheroid at zero temperature. 
As explained above, the penetration depth (5goid in gold of 
the electromagnetic fields that contribute to the potential 
U is always less than 20nm, for separations d larger than 
Xp and less than a few microns. A nano-particle of char¬ 
acteristic size satisfying the condition ^gold ^ ^ ^ d, 
can be modeled as perfectly reflecting. For such a parti¬ 
cle, both the electric and magnetic dipolar polarizablities 
need to be considered. The electric dipolar polarizability 
is given by Eq. (12) with e ^ oo. The dipolar mag¬ 
netic polarizability cP^ coincides with that of perfectly 
diamagnetic spheroid, and can thus can be obtained by 
setting /i = 0 in the formula for the magnetic polarizabil¬ 
ity of a magnetizable spheroid, given by 


FIG. 7: Stability diagram for a gold oblate spheroid (“pan¬ 
cake”) with ns = 0.7 at T = 0 K. 


M ^ V 

"" 47r 1 + (m - l)n, ■ 


(15) 


in Fig. Only the x and y axes occur as stable direc¬ 
tions, and the boundaries of the stable regions are simply 
given by Ri = R 2 and i7i ^ oc. A “needle” at T = 300 
K is always oriented along the z-axis in the parameter 


Since the dipolar polarizabilities of a perfectly con¬ 
ducting particle are frequency independent, the fre¬ 
quency integrals in Eq. @ can be performed analytically. 
At T = 0, the potential is then given by the explicit ex¬ 
pression 




































u = - 


hcv 


327r2n3 (1 — nl) 

11 + 693n3 - 266n| (f 


1 + 9n3 - 


1 - 3n3 


210 

25 + 38n3 dP 
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- 1 , 


r cos(26>) + f ^-^ 

^ 1^2 / V ^2 



This result shows again clearly that for a flat surface there 
is no orientation dependence of the potential. In Fig. 
we show the stable orientations of a prolate spheroid with 
ns = 0.7 (“pancake”), as a function of d/Ri and 
In Fig.j^an analogous plot is shown for a prolate spheorid 
with ns = 0.2 (“needle”). It is interesting to note that 
there exist values of RxjR^ for which the stable orienta¬ 
tion changes as the distance d is varied. For a spherical 
surface {Ri = R 2 ) ^ “pancake” prefers to sit parallel to 
the surface (symmetry axis oriented along the z axis) if it 
is located inside the sphere (negative radii of curvature) 
while a “needle” points towards a spherical surface if is 
outside the surface (positive radii of curvature). 



d/Ri 


FIG. 9: Stability diagram for a gold oblate spheroid (“pan¬ 
cake”) with n 3 = 0.7 and Ri = 20 //m at T = 300 K. 


Finite temperatures modify the stability diagrams: We 
consider again room temperature, T = 300 K, and as¬ 
sume that Ri = 20 /im (which we have to specify here 
explicitly since it is compared to the thermal wave length, 
contrary to the T = 0 case). As can be observed from 
Figs. thermal fluctuations reduce the stability re¬ 

gion for orientations of a “pancake” along the z-axis while 
increasing the stability for z-axis orientations of a “nee¬ 
dle.” In the latter case there is a change of the pre- 



FIG. 10: Stability diagram for a gold prolate spheroid (“nee¬ 
dle”) with n 3 = 0.2 and Ri = 20 /im at T = 300 K. 


ferred orientation for almost all ratios R 1 /R 2 with in¬ 
creasing distance d from either x- or //-orientation to a 
^^-orientation. 


IV. CONCLUSIONS ^ OUTLOOK 

On symmetry grounds it is expected that the Casimir- 
Polder force on an anisotropic particle, characterized 
by electric and magnetic dipolar polarizability tensors, 
should depend on its orientation relative to a nearby sur¬ 
face, with a torque rotating the object to energetically 
favorable alignment. Actually, for perfect conductors at 
zero temperature, and asymptotically at large distances, 
the interaction depends only on the trace of the static 
polarizability tensor, and orientation dependence at large 
separations is generically weak. At short distances, com¬ 
parable to the size of the object, strong orientation de¬ 
pendence is inevitable, selecting a favorable alignment 
for contact. For example, a prolate spheroid (pancake) 
will position itself with symmetry axis perpendicular to 
a flat surface {z direction), while an oblate (cigar) one 
































9 


will have its axis parallel to the surface {{x^y) plane). A 
curved surface, with distinct radii of curvature, will then 
break the rotational degeneracy of the oblate spheroid 
parallel to the plate. 

In this paper, we have studied the effects of surface cur¬ 
vature for the Casimir-Polder force on anisotropic nano¬ 
particles. The gradient expansion holds in an intermedi¬ 
ate range of separations, larger than the particle size, but 
smaller than the radii of curvature. While the expressions 
we find are quite generally valid- for arbitrary polariz¬ 
ability tensors and general material properties- we have 
focused on the easily visualizable case of spheroids near 
gently curved perfect conductors. We find that the inter¬ 
play of surface curvature and particle anisotropy leads to 
an orientation dependent interaction which is quite sensi¬ 
tive to temperature, separation, and dielectric response. 
While the minimum energy orientation is either perpen¬ 
dicular to the surface, or aligned to one of principal axes 
of curvature, the preferred alignment can change with 
temperature or separation to the surface. 


It should be noted that the computed orientation- 
dependence is a small fraction of the net Casimir-Polder 
interaction, complicating potential experimental probes: 
Freely suspended particles will be absorbed by the sub¬ 
strate, while trapped particles need to be cooled to very 
low temperatures before orientation preferences can be 
manifested. Nevertheless, it has been suggested [29] that 
such forces may be implicated in absorption properties 
of anisotropic molecules. A quantum treatment of the 
problem, applicable to the scales of molecular adsorp¬ 
tion, would thus be a valuable extension. 


Acknowledgments 

We thank R. L. Jaffe for valuable discussions. This 
research was supported by the NSF through grant No. 
DMR-12-06323. 


[1] T. A. Pasquini, M. Saba, G.-B. Jo, Y. Shin, W. Ketterle, 
D. E. Pritchard, T. A. Savas, and N. Mulders, Phys. Rev. 
Lett. 97 , 093201 (2006). 

[2] H. Oberst, D. Kouznetsov, K. Shimizu, J. I. Fujita, and 
F. Shimizu, Phys. Rev. Lett. 94, 013203 (2005). 

[3] B. S. Zhao, S. A. Schulz, S. A. Meek, G. Meijer, and W. 
Scho llkopf, Phys. Rev. A 78 , 010902(R) (2008). 

[4] J. D. Perreault, A. D. Gronin, and T. A. Savas, Phys. 
Rev. A 71 , 053612 (2005); V. P. A. Lonij, W. F. Holm¬ 
gren, and A. D. Gronin, ibid. 80 , 062904 (2009). 

[5] J. B. Pendry, J. Phys.: Gondens. Matter 9, 10301 (1997). 

[6] S. A. Biehs and J. J. Greffet, Phys. Rev. B 81, 245414 

( 2010 ). 

[7] O. A. Vasilyev, E. Eisenriegler, and S. Dietrich, Phys. 
Rev. E 88, 012137 (2013). 

[8] S. Kondrat, L. Harnau, and S. Dietrich, J. Ghem. Phys. 
131 , 204902 (2009). 

[9] H.B.G. Gasimir and D. Polder, Phys. Rev. 73 , 360 
(1948). 

[10] G.L. Klimchitskaya, U. Mohideen, and V.M. Mostepa- 
nenko. Rev. Mod. Phys. 81 , 1827 (2009). 

[11] Casimir Physics^ edited by D.A.R. Dalvit et al, Lecture 
Notes in Physics Vol. 834 (Springer, New York, 2011). 

[12] B. V. Derjaguin and LI. Abrikosova, Sov. Phys. JETP 3 , 
819 (1957); B. V. Derjaguin, Sci. Am. 203 , 47 (1960). 

[13] K. A. Milton, E. K. Abalo, P. Parashar, et al., Phys. Rev. 
A 91 , 042510 (2015). 

[14] T. Emig, N. Graham, R. L. Jaffe, and M. Kardar, Phys. 
Rev. A 79 , 054901 (2009). 

[15] R. Balian and B. Duplantier, Ann. Phys. (N.Y.) 104 , 300 
(1977); 112, 165 (1978). 


[16] A. Lambrecht, P. A. Maia Neto, and S. Reynaud, New 
J. Phys. 8, 243 (2006). 

[17] T. Emig, N. Graham, R. L. Jaffe, and M. Kardar, Phys. 
Rev. Lett. 99 , 170403 (2007). 

[18] R. Messina, D.A.R. Dalvit, P. A. Maia Neto, A. Lam¬ 
brecht, and S. Reynaud, Phys. Rev. A 80 , 022119 (2009). 

[19] V.B. Bezerra, E.R. Bezerra de Mello, G.L. Klimchitskaya, 
V.M. Mostepanenko, and A.A. Saharian, Eur. Phys. J. G 
71 , 1614 (2011). 

[20] G. D. Eosco, E. G. Lombardo, and E. D. Mazzitelli, Phys. 
Rev.D 84, 105031 (2011). 

[21] G. Bimonte, T. Emig, R. L. Jaffe, and M. Kardar, EPL 
97 , 50001 (2012). 

[22] G. Bimonte, T. Emig, and M. Kardar, Appl. Phys. Lett. 
100 , 074110 (2012). 

[23] V.A. Golyk, M. Kruger, A.P. McGauley, and M. Kardar, 
EPL 101 , 34002 (2013). 

[24] G. D. Eosco, E. G. Lombardo, and E. D. Mazzitelli, Phys. 
Rev. A 88, 062501 (2013).. 

[25] We normalize the waves |k, Q) as in Ref. [27]. Note 
though that the choice of normalization is irrelevant for 
the purpose of evaluating the trace in Eq. Q. 

[26] G. Bimonte, T. Emig, and M. Kardar, Phys. Rev. D 90 , 
081702(R) (2014). 

[27] A. Voronovich, Waves Rand. Media, 4 , 337 (1994). 

[28] D. B. Hough, and L. R. White, Adv. Golloid. Interface 
Sci. 14 , 3 (1908). 

[29] P. Thiyam, P. Parashar, K. V. Shajesh, G. Persson, M. 
Schaden, 1. Brevik, D. E. Parsons, K. A. Milton, O. 1. 
Malyi, and M. Bostrm, arXiv: 1506.01673 (2015). 



